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Abstract. The introduction of a steady-state spatially-periodic Raman coupling 
between two components of an ultracold atomic gas produces a dressed-state gas 
with an anisotropic and tunable dispersion relation. A Bose-Einstein condensate 
formed in such a gas is consequently characterized by an anisotropic superfluid critical 
velocity. The anisotropic dissipation of superfluid flow is quantified by considering 
the scattering of impurities flowing through this superfluid. A gradual transition 
from the isotropic nature of an uncoupled Bose-Einstein condensate to the anisotropic 
periodically-dressed condensate is obtained as the strength of the Raman coupling is 
varied. These results present a clear signature for future experiemental realizations of 
this novel superfluid. 
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The experimental attainment of quantum degenerate dilute gases, composed of 
both bosonic and fermionic atoms, has created many new opportunities in the study of 
quantum fluids. In particular, recent years have seen a flurry of experimental studies of 
superfluidity in dilute, scalar Bose-Einstein condensates. These include the observations 
of critical velocities for superfluid flow about microscopic |J and macroscopic [[J 
obstacles, the onset of turbulent flow above this critical velocity 0], beautiful studies of 
quantized vortices and vortex lattices in rotating Bose-Einstein condensates f|, || P , and 
other manifestations of irrotational superfluid flow |7], || . These experiments are closely 
analogous to those performed on another scalar superfluid, liquid 4 He, but apply new 
experimental probes and allow new insights due to the vastly different parameter regime 
of the dilute atomic gases, and their amenability to a new set of tools for manipulation 
and probing. 

Dilute atomic gases also offer the opportunity to create novel quantum fluids by 
using these various tools to manipulate the internal and external states of the atoms 
comprising the superfluid. For example, different types of mult i- component condensates 
have been studied: externally-coupled two-component condensates of 87 Rb at JILA 
(reviewed in ||), and F = 1 spinor condensates of sodium at MIT (reviewed in ||10|| ). 



One experiment at JILA explored how a spatially-selective coupling between two trapped 
hyperfine levels of a rubidium condensate leads to a "winding" of the order parameter, 
akin to the phase winding of the order parameter which occurs due to rotation about a 
vortex core, and how this winding is recurrently "undone" due to the motional dynamics 
of a two-component condensate O] . This observation illustrates that vortices in such a 
novel system should not be metastable (similarly predicted for spinor condensates ||12||), 
a major modification from the superfluid behaviour of a single component gas. 

In a similar vein, we have previously considered the novel superfluid properties 
introduced to a two-component Bose condensed atomic gas which is coupled by a 
spatially periodic coupling field This spatially-periodic coupling results from 

intersecting, non-collinear laser beams which induce a Raman coupling between two 
internal states of the ultracold gas. The Raman coupling introduces an anisotropy to this 
periodically-dressed fluid which should manifest itself in an anisotropic superfluid critical 
velocity. Controlling various parameters of the Raman coupling laser beams allows one 
to dramatically alter this critical velocity, leading to the dissipation of superfluid flow 
at velocities much lower than the critical velocity for the uncoupled Bose condensate. 

In this paper, we develop further a theory describing the dissipation of superfluid 
flow in a periodically-dressed Bose-Einstein condensate. Using the Bogolibov 



approximation theory which was introduced in Ref. |13| , the possible disspation of 
superfluid flow is treated by analyzing the scattering of massive impurities which flow 
through the superfluid. In analogy with similar calculations for a scalar Bose-Einstein 
condensate jl^, [I], impurity scattering at low impurity velocities is found to be 
suppressed by the energetics of the quasi-particle spectrum (related to the Landau 
criterion for superfluidity [PSfl), by the suppression of density fluctuations in the low- 
momentum phonon regime |17| , and by variations of the internal-state compositions 
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of quasiparticles at different wavevectors. These effects combine to give a smooth 
evolution in the disspation rates for superfluid flow as the Raman coupling strength is 
increased, gradually converting an uncoupled two component Bose-Einstein condensate 
(essentially a scalar superfluid) to a periodically-dressed condensate with dramatically 
different superfluid properties. This superfluid dissipation rate gives an appealingly 
clear signature for future experimental studies. 



1. Quasi-particle modes of the periodically-dressed condensate 



As presented in Ref. ]T3), a periodically-dressed Bose-Einstein condensate (PDBEC) 
is formed of atoms of mass m with two internal states, A and B, separated by an 
energy difference hujQ. These atoms are exposed continuously to two laser beams of 
wavevectors ki and k 2 and frequencies lo\ and u 2 , which couple the states A and B via 
a Raman transition (see Figure In such a Raman transition, an atom in internal 
state A with wavevector q — k/2 absorbs a photon from beam 1 and emits a photon 
into beam 2, arriving in internal state B with wavevector q + k/2 gaining a momentum 
hh = h(hi — k 2 ) and a kinetic energy K6 = h{uj\ — u> 2 — u; ). 



(a) A JiV |e> < b > 




k — k,| - k 2 



k 2 , ct>2 



Figure 1. Engineering propeties of a periodically-dressed atomic gas. (a) Laser beams 
of frequency oj\ and u>2 induce Raman transitions between internal states \a) and \b). 
The beams share a common large detuning A from the excited atomic state, and this 
state is adiabatically eliminated in the theoretical treatment. The lasers are detuned 
by an amount 5 from the two-photon Raman transition, (b) Such a Raman transition 
imparts a momentum transfer of Kk = h(ki — k^), where ki and k2 are the wavevectors 
of the Raman coupling lasers. 

Keeping in mind this coupling between atoms of different momenta, we let a q (a q ) 
be the creation (annihilation) operator for an atom in internal state A and wavevector 
q — k/2, and similarly let 6 q (6 q ) be the creation (annihilation) operator for an atom in 
internal state B and wavevector q + k/2. We refer to these as the bare-state creation 
and annihilation operators. The justification for this particular notation lies in the 
dressed-state picture |18[], in which one considers the driving optical fields to be part 
of the quantum system. In this picture, the Raman transition connects the states 
|a q ) = \A, q— k/2; N± + 1, N 2 ) and |6 q = \B, q + k/2; N%, N 2 + 1), where |a q ) represents 
an atom in state A and wavevector q — k/2 and an optical field with Ni + 1 photons in 
beam 1 and N 2 photons in beam 2, and similarly for |6 q ). In the dressed-state picture, 
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|a q ) and |6 q ) are states for which the total momentum of the atom + photon system 
is equal - by a proper choice of intertial frame, this total momentum is taken to be 
hq). The states are separated by an energy difference hS, and coupled by a static and 
(in the dressed-state picture) spatially uniform Raman coupling. The system is then 
represented by a many-body Hamiltonian of the form 

,2 / , \ 2 




Here Q is the two-photon Rabi frequency characterizing the strength of the Raman 
transition, determined by the polarizations and intensities of the driving laser beams 
and by dipole matrix elements of the atoms. 

This Hamiltonian is diagonalized by transforming to the dressed-state creation and 
annihilation operators /i q , /x q , 7r q , 7r q which are defined as linear combinations of the 
bare-state operators 

[b^ J { -sin^/2 cos0,/2 / V ) 

The mixing angle 9 q is defined by the relation tan# g = to/ {8 — h<\ ■ k/m). The dressed 
states have energies 



(2) 




Here, the minus sign refers to the lower dressed-state dispersion curve and to the 
operators fi n and /A, while the plus sign refers to the upper dressed-state dispersion 
curve and to the operators 7r q and 7r q . 

A Bose-Einstein condensate formed from this periodically-dressed atomic gas 
has a macroscopic population of N Q atoms in the lowest energy momentum state: 
this state lies on the lower dressed-state dispersion curve and is taken to have 
momentum HQ, |19|]. In Ref. ||13|| , a Bogoliubov-approximation theory was developed 



to account for the effects of weak interatomic interactions on a PDBEC. An interaction 
Hamiltonian of the form 7Y m ^ = | X) q ( n q n -q — N) was considered, where N is the 
total number of atoms in the system and n q = J2k ( a k+ q a k + ^k+ q ^k) is the Fourier 
transform of the density operator. Quasi-particle energies and operators were found 
by diagonalizing a 4 x 4 matrix (-ffq)ij = {£q)ij + ^Xi%j where £ q is a diagonal matrix 
with entries (£~ , — £r q ,£ q , — £t q ) with £ q = ^(u;Q +q — uJq), \i = gN is the chemical 
potential, and x = (cos A q , — i cos A_ q , sin A q , — ism A_ q ), where A q = (#Q +q — 0q)/2. 
Diagonalization yields 

tf q = M q # q M q 1 (4) 

where Hq is a diagonal matrix with entries fr(d)Q +q , — u>Q_ q , WQ +q , — ^Q- q ) by which the 
lower (minus sign) and upper (plus sign) quasi-particle energies Uu)q + are defined. One 
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also obtains the creation and annihilation operators for the lower (/z q , /z q ) and upper 



7T n , 7T 



branch quasi-particles through the following relation: 



/ a Q+q \ 

t 



%a, 
be 



Q+q 



/ /^Q+q \ 

t 



# q M q 



/ /^Q+q 
^Q-q 



7TQ+q 
V ^Q-q 



(5) 



^Q-q 
^Q+q 

Here _R q contains trigonometric functions of #Q +q /2 and #Q_ q /2 in accordance with Eq. 
0. The Bose commutation relations take the form 

[ i^j 

Xi, Xj\ = Aj = S 1 i = j = 1 or i = j = 3 (6) 
1—1 i = j = 2 or i = j = 4 

where x represents any of the four- component vectors in Eq. |5|. These relations are 
clearly preserved by the rotation matrix, i.e. R^AR^ = A, while the imposition of Bose 
commutation relations for the quasi-particle operators enforces a normalization criterion 
M q AMt = A. 

Experimentally, values for the Raman detuning 5, the Rabi frequency Q, the 
chemical potential [i and the Raman momentum transfer hk can be chosen over a wide 
range of parameters. For example, let us consider a particular realization of a PDBEC 
using two hyperfme states of 87 Rb. One may consider the case of Raman excitation 
using two lasers both tuned near the D2 optical transition at a wavelength A ~ 780 
nm. Thus, the Raman recoil energy can be readily tuned in the range 

< Ek < h x 15 kHz. Properties of a PDBEC are governed by the ratios hS/E^, Vt/E^ 
and fi/Ek- The Raman detuning 5 and the Rabi frequency Q can be chosen nearly 
arbitrarily, reaching maximum values which are orders of magnitude larger than E^. 
The chemical potential jx = Anh 2 an j m has a typical value of 1 kHz for a density of 
n~lx 10 14 cm -3 and a scattering length a ~ 100 bohr - thus the ratio fi/Ef. can be as 
small as ~ 1/10 for counter-propagating Raman excitation lasers, and greatly increased 
by using Raman beams intersecting at small angles. 

For illustration, we show in Figure |2| the quasiparticle dispersion relations for a 
PDBEC with a chemical potential [i/Ej. = 1, a Raman detuning of of Tib~/Ek = 1/2 
and Raman coupling strength hQ/Ek varying between and 2. For this case of 5 > 0, 
a condensate is formed in the Q ~ — k/2 state of the lower dispersion relation. In 
the absence of Raman coupling, this equality is exact, and the condensate is formed in 
the \B) internal state. Quasiparticle excitations in the \B) state follow the standard 
Bogoliubov result, while quasiparticles in state \A) have a free-particle dispersion 
relation - this can be though of as due to the lack of a spatial interference term between 
the condensate and excitations in the \A) state at any wavevector. For small Rabi 
coupling strengths, one speaks more properly of excitations in the lower and the upper 
quasiparticle dispersion relations, although the energies and internal-state compositions 
of these states remain relatively unchanged for wavevectors q far from the Doppler- 
shifted Raman resonance condition 5 — q • k/m = 0. The lower quasiparticle dispersion 
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relation has a local minimum near q = k/2 caused by the minimum of the dispersion 
relation for the bare \A) state atoms. By analogy to the roton minimum in the excitation 
spectrum of superfluid 4 He JT^] , we denote this feature as an "artificial roton," although 
we mean to suggest no similarity between the specific structure of a roton and of the 
plane-wave excitations of a PDBEC. As the Rabi frequency is increased further, a 
greater mixing of the internal state composition of the Bose condensate and of the 
quasiparticle excitations occurs, and energy levels are shifted. In particular, the energy 
of the "artificial roton" is shifted upwards due to interactions as the internal state 
composition of the Bose condensate and of these excitations becomes more similar. 




Figure 2. Quasiparticle dispersion relations for a PDBEC. Shown are the energies of 
quasiparticle in the lower (created by £tQ+ q ) and upper (created by 7Tq +q ) dispersion 
curves with momentum q with respect to the condensate. Excitations parallel to the 
Raman momentum transfer are considered, with the definition q = q • k/fc. Black 
curves represent PDBEC 's with HS/Ek = 1/2, fj, = 1 and Ml/Ek — 1/2 (dashed), 1 
(dotted) and 2 (solid). The gray dot-dashed curve shows the Bogoliubov excitation 
spectrum for a single component condensate with (i = 1. 



2. Impurity scattering 

We now consider the scattering of impurities passing through a uniform PDBEC. 
We consider impurities of mass M c which have a free-particle dispersion relation 
E c = h 2 k 2 /2M c where k c is the impurity wavevector. Such impurities interact with 
the two components of the periodically-dressed condensate through an interaction of 
the form 

T^SCat — ^ X! ( a L+qflmCi +q Cl + &L+q&mCi+q C l) (7) 
m,l,q 

Here a q , a q , 6 q , and 6 q are defined as above, and c q and c q are the creation and 
annihilation operators, respectively, for the impurity particles. The coupling strength 
for the impurity scattering is A, which we parameterize as A = 2ith 2 a c /M' . This 
parameterization correctly describes low-energy s-wave scattering with a scattering 
length of a c and the reduced mass M' = mM/ (m + M) . Note that we have made 
the assumption that the scattering strengths between the impurity and atoms in states 



Dissipation of superfluid flow in a periodically-dressed BEC 



7 



A and B are equal, and assume further that state changing collisions do not occur (the 
impurity is considered "non-magnetic"). We may consider also the limit as the impurity 
mass becomes infinite while the impurity velocity v = hk/m remains constant (k is the 
impurity wavevector). We thereby obtain a description of the PDBEC flowing past a 
rigid obstacle with a velocity v, and the concept of impurity scattering is replaced by 
the dissipation of superfluid flow. 

The scattering rate of these impurities can be calculated using a perturbative 
approach as follows [|I^,[l4]]. We consider the system to be initially in the state |0)As|k)c, 
i.e. a PDBEC with no excitations, and a single impurity particle of wavevector k. The 
scattering Hamiltonian of Eq. [7] may couple this state to the final state |/)AB|k — q)c 
in which the periodically-dressed gas is in the state |/) and the impurity wavevector 
changes from k to k — q. The matrix element for this transition is A(/|n q |0), where 
n q is described above, and where the AB subscript labeling the quantum states of 
the periodically-dressed gas has been dropped. We then obtain the rate of impurity 
scattering by summing over all final states |/) and momentum transfers q using Fermi's 
Golden Rule. 

We now make use of the Bogoliubov approximation. If we consider only weak 
scattering, we may restrict the final states to those containing a single quasiparticle, i.e. 
\f) = Iq"* 1 ) where |q~) = /iQ +q |0) and |q + ) = 7rQ +q |0). The rate of impurity scattering 
is then 

1 = fr A2 E S M) 6 (^Q+q - v • q + ^-j + S n {q)5 - v • q + ^ 

where v = hk/M is the impurity velocity. We have introduced the quantities S^q) and 
S n (q) which relate to the static structure factor describing density-density correlations 
in the Bose condensate. For a single-component Bose-Einstein condensate, the static 
structure factor S(q) has been evaluated and measured |2(| ^T[. For the PDBEC, 
we separate the contributions of the lower and upper excitation spectra as 

^(q) = ^(0|n_ q |q-)(q> q |0> (9) 
iV 

^(q) = ^(0|n_ q |q + >(q + |n q |0> (10) 

iV 

Physically, these quantities describe two effects. First, the magnitude of S M and S w 
describes the contributions of quasiparticles in the lower (/i) or upper (tt) quasiparticle 
state to density fluctuations (as opposed to phase fluctuations, which dominate, for 
example, for phonon excitations of a scalar Bose condensate pO" , 21 , |23|). Second, these 
quantities describe the degree to which the internal state composition of the condensate 
matches that of the quasiparticle. For instance, a condensate which is predominantly 
in the \B) internal state will not scatter strongly into quasiparticle states in which the 
\A) internal state is dominant, and thus the structure factor describing this scattering 
will be small. 

Before proceeding to a calculation of and S n , it is useful to consider the scattering 
rate Tq for an ideal-gas, single-component Bose-Einstein condensate as the target gas. In 
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this case, we may conveniently take S^q) = 1, ^(q) = and u = hq 2 /2m. Replacing 
J2d = V/i^rr) 3 f d 3 q, and substituting for A we evaluate Eq. || and obtain simply 



X 2 N VM ,2 v N 2 

T = — 4 = -x 4vra c x v, 

Ten V 



(11) 



i.e. the conventional expression for the scattering rate where N /V is the density of the 
target gas and Aira^ is the collision cross section. We may then determine the impurity 
scattering rate in a PDBEC relative to that in an ideal gas, i^ot = r/To = F fJt + F n : 

h 2 hq 2 



4vr,M^/ d3q ^ (q)5 r Q+q ' 2M C 

h 2 r .« „ , v „ / _ ^ hq 2 



- v • q 



2M, 



- vq 



(12) 
(13) 



c 



3. Structure factor formalism 



To calculate the structure factors and S n we start by expressing the density operator 
n q in terms of the quasiparticle operators 



a 



t 

Q+m+q 



-2GlQ-m-q ^Q +m+q 



-ib 



Q m — q 



t 

Q m 



HI 



o ( Acj+m+q ^Q-m-q 



^Q+m+q — *^Q-m-q ) ^q+m-Rq+m-Rm^m 



(14) 



Acting on a PDBEC, the only important terms in the sum are those which include the 
condensate operators /}q and /2q. Applying the relations jiQ +q \0) = 7TQ+ q |0) = for 
q 7^ 0, and the relation M = 1, we define matrices a = M^R^Ro and (3 = i?o-R- q M_ q 
and obtain 

Wo 



^ql ) = - ^ (("11 + ^12 + *#L2 + ^22)^10) + (CK31 + Za 32 + «/3l4 + /3 24 )vT q |0)) (16) 



So 



1 2 
•^(q) = ^ Kl + ^12 + ^12 + Ad 

1 i2 

SV(q) = 1 |«3i + 2a 32 + ifiu + #24 1 



;i7) 



The static structure factors S^(q) and S n (q) for a PDBEC are shown in Figures 
|3] and The specific values for the Raman detuning 5, the Rabi frequency Q and the 
chemical potential fi are relevant for a particular realization of a PDBEC using two 
hyperfine states of a 87 Rb, as discussed above. The division of the scattering weight 
between the lower (S^) and upper (S n ) quasiparticle branches can be understood by 
considering the composition of the dressed-state quasiparticles, which is related by 
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Eq. to the mixing angle # q . This angle ranges from # q = for q ■ k — > — oo to 
# q = 7i for q • k — ► +oo, reaching # q = tt/2 at the Doppler-shifted Raman resonance 
5 — q • k/m = 0. For small Rabi frequency, # q , and hence the internal state composition 
of the dressed states varies rapidly about the Raman resonance condition, while for 
large Rabi frequency, this variation occurs over a broader range of q. This behaviour 
is reflected in the structure factors, as illustrated by Figure [| By continuity, the 
internal state composition of the condensate in the hq state is nearly identical to that of 
excitations in the lower quasiparticle branch (created by p.Q +q ), and nearly orthogonal 
to that of excitations in the upper quasiparticle branch (created by vrQ +q ). Thus, since 
elastic scattering from a "non-magnetic" impurity does not affect the internal state 
distribution, S n (q) — generally for small q. Away from q = 0, the division of the 
scattering weight between the upper and lower quasiparticle branches depends on the 
strength of the Rabi coupling. For small Rabi frequencies, the scattering strength shifts 
abruptly from the lower to the upper quasiparticle branch at the Doppler-shifted Raman 
resonance. For large Rabi frequency, the changes in the scattering strengths are more 
gradual, with scattering into both branches showing a significant weight for a wide range 
about the Raman resonance condition. This variation is particularly relevant for the 
scattering of impurities near the Landau critical velocity, as discussed below. 




Figure 3. Variations of the static structure factors (a) S^q) and (b) S 7r (q) with the 
Raman coupling strength The wavevector q is colinear with the Raman momentum 
transfer k with q = q- k. Three black curves describe a PDBEC with hS/Ek = 1/2, 
fx/Ek = 1, and hfl/E k — (solid), 1/4 (dashed) and 1 (dotted). Scattering for q ~ 
and q < occurs only into the lower quasiparticle branch and for q 3> k into the upper 
quasiparticle branch. In other regions, stronger Raman coupling leads to eigenstates 
with mixed populations in the internal states \A) and \B) over a broader range of q, and 
thus impurity scattering occurs into both the lower and upper quasiparticle branches. 
Also shown is the structure factor 5(q) for a scalar Bose-Einstein condensate (gray 
dot-dashed curve) with fx/Ek = 1. 

As illustrated in Figure f|, interatomic interactions suppress scattering at small 
momentum transfers. For a single-component homogeneous Bose-Einstein condensate, 
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Figure 4. Static structure factor (a) 5 M (q) and (b) S 7r (q) for a PDBEC. The 
wavevector q is colinear with the Raman momentum transfer k with q — q ■ k. Three 
black curves describe a PDBEC with %5/Ek = 1/2, Ml/Ek — 1, and fi/Ek = (solid), 
1/4 (dashed) and 1 (dotted). The condensate in the hq state is predominantly in 
the internal state B; thus scattering predominantly occurs into the lower quasiparticle 
branch for q < and into the upper quasi-particle branch for large q > 0. For q ~ 0.5 
both S fl (q) and S- K (q L ) are significant since quasiparticles contain large fractions of 
both | A) and \B). Interactions suppress scattering at small momentum transfer, an 
effect seen also in the structure factor S(q) for a scalar Bose-Einstein condensate (gray 
dot-dashed curve) with (x/E). = 1. 



the structure factor S(q) = q / Jq 2 (q 2 + 2£ 2 ) is small for wavevectors of magnitude 
q < l/£ where £ = i/2m/i//j 2 is the healing length. Similarly, for the PDBEC, we find 



a suppression of the structure factors S^q) and S 7T (q) for q < kJ2[i/Ek, as shown 
in the figure. This suppression of the structure factor in a PDBEC could be probed 
experimentally as it was for a scalar Bose-Einstein condensate by studying shallow- 



angle optical Bragg scattering [2C. 21 



4. Dissipation of superfluid flow 

Finally, let us use the formalism developed above to describe superfluid properties 
of a PDBEC. In Ref. ||13||, we applied the Landau criterion to determine the critical 



velocity for superfluidity in a PDBEC given the quasiparticle dispersion relations. 
This criterion gives the superfluid critical velocity along a direction q as having a 
magnitude v c = min 9]i Ei(qq)/hq where E^q) are the energies of possible excitations 
with wavevector q; using the momentum definitions for a PDBEC, these would be the 
quasiparticle energies ft,cjQ +q . In an unperturbed, single component weakly-interacting 
condensate, this critical velocity t>° is equal to the Bogoliubov speed of sound c s = \j\ijra 
in all directions. In a PDBEC, the anisotropic quasiparticle dispersion relation leads 
directly to the prediction of an anisotropic superfluid critical velocity, and to the 
possibility of a critical velocity lower than that of the unperturbed condensate. 
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One may determine this critical velocity graphically by finding the line of minimum 
absolute slope which connects the origin and the quasiparticle dispersion curve (see 
Figure ^). For the PDBEC, the superfluid critical velocity may be suppressed by the 
presence of the "artificial roton" feature in the lower dispersion curve. This suppression 
occurs in the direction of the Raman momentum transfer Kk for the case of 5 > (as 
shown in the figure), or in the opposite direction for 5 < 0. 

However, while the Landau criterion determines the minimum velocity for the 
dissipation of superfluid flow, it does not quantify the onset of such dissipation as the 
velocity is increased beyond the critical value. Indeed, for extremely weak Raman 
coupling, even though the superfluid critical velocity is strictly different, one should 
expect the fluid properties of a PDBEC to be nearly identical to that of a two-component 
BEC in the complete absence of Raman coupling. We can quantitatively assess this 
description using the formalism developed above for the determination of the matrix 
elements and S^, and therefrom the normalized scattering rates F^ and F n . 

The results of such calculations for infinitely massive impurities (i.e. superfluid 
flow about rigid obstacles) are presented in Figures | and || Values of and F n 
were determined by numerical integration in which the delta function in Eqs. O and 



13 was approximated by a narrow Gaussian distribution whose width was adjusted to 



give reasonable numerical convergence. As we expect, for small Rabi frequencies Q the 
properties of a periodically-dressed BEC are quite similar to those of an uncoupled two 
component Bose-Einstein condensate. The scattering probability for impurity atoms 
(related to the viscosity of the fluid) is indeed non-zero at velocities lower than the 
superfluid critical velocity for the unperturbed gas, demonstrating that the critical 
velocity is suppressed along the direction of Raman scattering. This scattering at low 
velocities occurs through the excitation of quasiparticles near the local minimum (the 
"artificial roton") in the lower dispersion relation. However, this scattering channel is 
weak for small values of the Rabi frequency Q, and thus the scattering strength i^ot 
remains weak until the curve joins that of the uncoupled condensate and the scattering 
strength is increased by phonon scattering. 

The effects of increasing the Rabi frequency are twofold. First, as Q increases, the 
internal state composition of the "artificial roton" excitation becomes more similar to 
the internal state composition of the Bose-Einstein condensate. This causes an increase 
in the scattering strength (as indicated also in the calculations of shown in Figures 
f| and |3|). Simultaneously, because of this overlap of the internal state composition, 
the energy of these excitations is increased, approaching the energy given by the one- 
component Bogoliubov spectrum, albeit with an appropriately modified effective mass. 
Thus, the minimum velocity at which scattering occurs increases as the Rabi frequency 
increases. At moderate values of the Rabi frequency (where hQ/E^ ~ 1), the scattering 
rate vs. velocity has a two-step shape as scattering is first enhanced by strong scattering 
into the "artificial roton" state, and then later enhanced by phonon scattering when 
the velocity is near v®. As shown in Figure for one typical parameter setting, only F t 



contributes to the weak scattering at velocities smaller than v®. At velocities greater 
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Figure 5. Dissipation strength of superfluid flow in periodically-dressed Bose-Einstein 
condensates. The scattering rate F tot for infinitely massive impurities traveling through 
the condensate, normalized to the scattering in an ideal gas Bose-condensate at similar 
densities, is plotted vs. the velocity of the impurity. Here we consider velocities parallel 
to the Raman momentum transfer, with v = v • k/fc. Black curves show results for 
PDBEC's with varying strengths of the Raman coupling; for all curves, 6/Ek = 1/2 
and /i/Ek = 1 while huj/Ek = 1/2 (dashed), 1 (solid) or 2 (dotted). The gray dot- 
dashed curve shows results for a one-component condensate with fi/Ek = 1. The 
superfluid critical velocity of the PDBEC in the direction of the Raman momentum 
transfer is suppressed with respect to an uncoupled condensate. For weakly Raman 
coupling, the dissipation at low velocities is weak. This dissipation grows stronger as 
the Raman coupling strength is increased. 

This two-step curve gives a clear signature which should be visible in experimental 
probes of superfluid flow. Indeed, experiments such as those performed by Chikkatur 
et al 0, in which atoms of another ground state hyperfine level were used as scattering 
impurities in a dilute atomic Bose condensate, could be carried out in a straightforward 
manner. It would be interesting to consider further how other aspects of superfluidity, 
such as the onset of turbulent flow |2| or the metastability of quantized vortices [§], [5|, are 
affected by the anisotropic nature of a periodically-dressed Bose-Einstein condensate. 

In conclusion, dissipation of superfluid motion in a periodically-dressed Bose- 
Einstein condensate above the superfluid critical velocity was treated. A formalism 
was developed for the calculation of the static structure factors and S n describing 
density-density correlations from quasiparticles in the lower (/z) and upper (jr) excitation 
branches of the periodically coupled gas. Following a perturbation theory approach, 
the scattering rate of massive impurities, or equivalently the damping rate of superfluid 
motion past microscopic obstacles, was calculated. Raman coupling of a two-component 
condensed gas leads to a suppression of the superfluid critical velocity. The impurity 
scattering rate evolves smoothly as the Raman coupling strength is increased. The the 
superfluid damping is weak for weak Raman coupling, and becomes more significant as 
the Raman coupling strength is increased. 
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Figure 6. Normalized rates for production of lower (F M , dashed line) and upper (F n , 
dotted line) quasiparticle excitations due to scattering off infinitely massive impurities 
moving at velocities v in a PDBEC. Velocities parallel to the Raman momentum 
transfer are considered, with v = v • k/fc. Black curves show results for a PDBEC 
with H8/Ek = 1/2, fi/Ek = 1, and hfl/Ek — 1, with the solid curve showing the total 
normalized scattering rate -Ftot- The gray dot-dashed curve shows the normalized 
scattering rate for a one-component condensate with ji/Ek = 1. For v > (for this 
case where 5 > 0), impurity scattering at low velocities is enhanced as compared to the 
one-component condensate. This scattering leads to the creation of excitations in the 
lower dispersion curve (F M > 0) around the "artificial roton" feature in the quasiparticle 
dispersion relations. At velocities near the critical velocity for the unperturbed 
condensate, v®, scattering is allowed into the upper dispersion relation F^ > 0, and 
the normalized scattering rate approaches that of a one-component condensate. For 
v < 0, energy and momentum conservation imply that only scattering into the lower 
quasiparticle excitations is allowed at the velocities considered. 
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